The case for testing MOND using LISA Pathfinder 
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We quantify the potential for testing Modified Newtonian Dynamics (MOND) with LISA 
Pathfinder (LPF), should a saddle point flyby be incorporated into the mission. We forecast the 
expected signal to noise ratio (SNR) for a variety of instrument noise models and trajectories past 
the saddle. For standard theoretical parameters the SNR reaches middle to high double figures even 
with modest assumptions about instrument performance and saddle approach. Obvious concerns, 
like systematics arising from LPF self- gravity, or the Newtonian background, are examined and 
shown not to be a problem. We also investigate the impact of a negative observational result upon 
the free-function determining the theory. We demonstrate that, if Newton's gravitational constant 
is constrained not be re-normalized by more than a few percent, only contrived MONDian free- 
functions would survive a negative result. There are exceptions, e.g. free- functions not asymptoting 
to 1 in the Newtonian limit, but rather diverging or asymptoting to zero (depending on their mother 
relativistic MONDian theory) . Finally, we scan the structure of all proposed relativistic MONDian 
theories, and classify them with regards to their non-relativistic limit, finding three broad cases 
(with a few sub-cases depending on the form of the free function). It is appears that only the 
Einstein- Aether formulation, and the sub-cases where the free-function does not asymptote to 1 in 
other theories, would survive a negative result without resorting to "designer" free-functions. 



I. INTRODUCTION 

Einstein's theory of General Relativity (GR) and the 
ACDM standard model are two cornerstones of modern 
cosmology. They posit that the gravitational effects of 
large scale structures in the universe (such as galaxies 
and clusters of galaxies) cannot be explained by lumi- 
nous, baryonic matter alone, but rather that an addi- 
tional cold (pressureless) and dark (non-luminous) mat- 
ter component is needed. However, in the absence of 
direct observational evidence for dark matter, it remains 
nothing but a useful calculational device. For as long as 
this is true, it is scientifically healthy to explore alter- 
native explanations for the anomalous gravitational dy- 
namics, namely by modifying the theory of gravity itself. 

Modified Newtonian Dynamics (or MOND [1|) is one 
such scheme valid in the non-relativistic regime. It was 
first proposed to explain observed dynamical properties 
of galaxies without invoking dark matter. More recently 
it has been incorporated into relativistic theories 0t9j 
following from the ground breaking proposal of "TeVeS" 
by Bekenstein [Tj]. Relativistic extensions are needed for 
reasons beyond logical completeness: they are required 
to explain, for example, phenomenology associated with 
lensing and cosmology (8j , where dark matter is also usu- 
ally employed. When the whole picture is assembled the 
conflict between MOND and dark matter leaves consid- 
erable scope for doubt over the interpretation of new as- 
trophysical and cosmological data. A fair comparison 
requires re-evaluating, within each approach, the whole 
set of assumptions underlying the new observations. For 
this reason the debate would benefit from a direct probe, 
in the form of a laboratory or Solar System experiment. 
This has been proposed in various forms, namely in plan- 
etary data (9| appealing to the exterior field effect [io[ ■ 

The fact that MOND predicts anomalously strong tidal 



stresses in the vicinity of saddle points of the Newtonian 
potential has been advocated as one such decisive direct 
test [Ilj. The forthcoming LISA Pathfinder mission [IJ] 
presents the perfect opportunity for its realization, as a 
preliminary feasibility study has demonstrated |13l Il4| . 
The purpose of this paper is to provide a detailed quanti- 
tative evaluation of the power of a MONDian saddle test 
using LPF, predicated upon a scenario where a mission 
extension is granted. The extension would involve redi- 
recting the spacecraft from Lag rang e point LI to a saddle 
of the Earth- Moon-Sun system [ij] once its nominal mis- 
sion at LI is completed. In establishing the scientific case 
our efforts in this paper are twofold. 

In the first part of the paper we propose some basic 
data analysis tools and evaluate their expected perfor- 
mance. These tools are an adaptation of the "noise- 
matched filters" employed in gravitational wave detec- 
tion 15|. Their implementation benefits from a major 
simplification: for a saddle test we do know the tem- 
plate's starting point in time. A number of pitfalls and 
potential systematics found in detection of gravitational 
waves are therefore expected to be absent. The filter's 
optimal signal to noise ratio (SNR) allows us to quan- 
tify with a single number the predicted outcome for any 
experiment. Assuming a "standard" MONDian theory, 
the unknowns reduce to the instrument performance (the 
noise properties) and the trajectory past the saddle (its 
impact parameter). For each of these we can condense 
the expected outcome of a LPF test in a single number: 
the forecast SNR assuming MOND is correct. 

Our central results are in Section IIIIl particularly in 
Figs, m and [3 where the optimal SNR is plotted against 
noise level and saddle impact parameter. In a cataclysmic 
scenario for instrument performance and saddle approach 
we'd still achieved SNR sa 5. For less pessimistic assump- 
tions, high double figures are easily reached. We examine 



the effect of the spacecraft speed as it flies past the saddle, 
showing that just about any typical speed will turn out 
to be optimal. This is due to a remarkable coincidence, 
spelled out in Section IIVI and in the concluding section 
of this paper. In Section IIVI we also show that possible 
systematic errors, such as self-gravity or the Newtonian 
background, are in fact harmless. 

In the second part of this paper, and complementing 
the work just described, we spell out the generality, or 
otherwise, of the conclusions in the first part, and ex- 
amine the implications of a negative observational result. 
Just how comprehensively would the failure to detect the 
predicted high SNRs rule out the MONDian paradigm as 
a whole? As explained in Section |TT] the large menagerie 
of proposed relativistic MONDian theories practically all 
reduce to the same non-relativistic limit as TeVeS, and 
virtually all theories fall into 3 categories. One then is 
left with a free-function, /i, and the question is, how much 
leeway does it provide for evading a negative result? In 
Section |V] we review previously proposed free functions, 
rewriting them under a unified notation. We then lay 
down conditions for what should be permissible simple 
free functions at the most basic level (by simple we mean 
with a minimum number of regimes and scales) . Briefly 
we require that: (I) The theory shouldn't renormalize 
Newton's gravitational constant by more than a few per- 
cent in the Newtonian regime; (II) The theory should 
predict the usual MONDian effects when the Newtonian 
acceleration drops below acceleration scale oq. These 
constraints were implemented in TeVeS and set the stan- 
dard for a viable theory with useful astronomical impli- 
cations. We show that all natural functions satisfying 
these conditions result in similar SNRs for an LPF sad- 
dle test, as long as the impact parameter is smaller than 
400 km. The only exceptions are free-functions with a 
divergence (and with the rest of their domain excised) 
such as those suggested in [T^ [I^l; however these may 
fall foul of existing Solar system constraints and other 
requirements |17H2(]j . We leave for a future publication 
a more detailed examination of these functions. 

Of course one may open the doors to free-functions 
with more structure: e.g. functions with three or more 
regimes instead of the minimal two. In Section I VII we 
quantify how contrived the free function fi would have 
to be, for the theory to survive a negative result. We 
find that, for bounded functions, only a /i turning from 
1 (Newtonian regime) into an intermediate power-law, 
/J, ex z", and only then into the MONDian /i ex z, would 
be viable. The intermediate n would have to be very 
different from 1 even with undemanding requirements on 
impact parameter and noise. Another possibility are the 
free- functions suggested in |2l| - [24| . for which galaxy ro- 
tation fits and Solar system constraints are combined to 
motivate more structured free- functions. All of these in- 
voke three regimes (and so have two implied scales) . Even 
then, in Section FVl Al we show that these would be within 
striking distance for a saddle test with LPF. Thus, with 
the exception of diverging fi, only "designer" /i would 



bypass a negative result. Although this conclusion is de- 
rived for TeVeS-like theories it might be more general. 



jtions and provisos, the 
and a particular case 
might be the only rel- 



With the usual honorable excei; 

Einstein-Aether formulation 

of Milgrom's bimetric theory 

ativistic realizations of MOND to realistically survive a 

negative result from a saddle test, a matter we'll examine 

in a future publication. 

We conclude that a LPF test has both the power to 
detect MOND with a high SNR should it be true, and to 
rule it out, should a negative result be obtained. 



II. MONDIAN THEORIES 

One can find in the literature a large number of rela- 
tivistic MONDian theories. It is important to note that 
their complexity and differences arise from the require- 
ment that they should explain relativistic phenomena 
(such as lensing and structure formation) without ap- 
pealing to dark matter. However, in the non-relativistic 
regime, almost all of them reduce to the non-relativistic 
limit of TeVeS, which will be the focus of this paper. 
There are exceptions, however, and we spell them out 
here. In general the large profusion of relativistic MON- 
Dian theories reduces to only 3 types of non-relativistic 
limits, which we'll label type I, II and III. 

• Type I In these theories the non-relativistic dy- 
namics results from the joint action of the usual 
Newtonian potential $jv (derived from the metric 
via goo ~ —(1 + 2$Ar)) and a "fifth force" field, 
(j), responsible for MONDian effects. The total po- 
tential acting on non-relativistic particles is their 
sum: 



$ = $ 



N 



(1) 



Whilst the Newtonian potential satisfies the usual 
Poisson equation: 



V^$ 



N 



AttGp 
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the field 4> is ruled by a non-linear Poisson equation: 



V • {n{z)V(l)) = kGp, 



(3) 



where, for convenience, we pick the argument of the 
free function function ^ as: 



(4) 



47r oo 



where k is a dimensionless constant and ao is the 
usual MOND acceleration. In Section |V] we will 
say more on admissible functions /i, but in general 
we require that /i — >■ 1 when z ^ 1 and /i ~ z for 
z <C 1. (We use letter z instead of y to prevent a 
common source of confusion in the literature; see 
Section ES) 



Type II In these theories we also have <J> = $Ar + 0, 
but now the field cj) is ruled by a driven linear Pois- 
son equation, whose source depends on the New- 
tonian potential. In order to facilitate comparison 
with Type I theories (as explained below) we write 
the equation for (f) in these theories as: 



47r 



V • {v{v)V<^n) 



(5) 



where the argument of free function v is given by 



^^ / K\2 |V$jv| 

V 47r / oo 



(6) 



and we require that i^ — > 1 when v ^ 1 and 
V ^ l/\/v for w <C 1. However, it is possible in some 
models that that ^ — >■ in the same limit, with 
qualitatively very different implications. To distin- 
guish these two cases we call the latter (z^ ^- 0) 
type II A theories and the former {v —>■ 1) type IIB 
theories. 

Type III This was the original non-relativistic 
MONDian proposal, derived from a non-relativistic 
action principle (the so-called AQUAL |25l]). Cru- 
cially, here non-relativistic particles are sensitive to 
a single field $ which satisfies a non-linear Poisson 
equation: 



V • (M(a;)V$) = AirGp 



(7) 



Again, /i is a free function with a suitably chosen 
argument: 



ao 



(8) 



so that /i — > 1 when a; ^ 1 and jl ^ x for a; <C 1. 

A couple of remarks on this classification are in order. 
Firstly, note that in some theories $ = "E^m + </>, where 
S « 1 is a factor evolving slowly on cosmological time 
scales. We have set S = 1 throughout, for simplicity. 
Secondly, it is obvious that the two equations ruling type 
II theories may be rewritten as a single equation, ruled by 
a redefined v. This is a cosmetic matter and is beside the 
point. The real matter is whether G is renormalized, with 
regards to the bare G and the cosmological one (a mat- 
ter we will presently discuss in detail). This is encoded 
in the limiting behaviour of the free function ruling the 
non-relativistic equations, whether they're written as a 
single field or two fields system. It has been argued that 
in some relativistic formulations of type II theories the 
bare G, that appearing in cosmology and the total G rul- 
ing the non-relativistic equation are the same. We call 
such theories type IIA, and for them z^ — ?> in the New- 
tonian limit. Otherwise we call then type IIB, with a 
G renormalization, and with z/ — ;> 1. As we will see the 
matter has crucial implications. 

Virtually all relativistic MONDian theories proposed 
in the literature fall into these categories. TeVeS, the 



pioneering relativistic MONDian theory Q, has type I 
limit, but Sanders' stratified theory 0] is also type I. Mil- 
grom's bimetric theory 0, Q can be either type I or type 
II, depending on details. Einstein- Aether theories 0, 0] 
are unique in that they have a non-relativistic limit of 
type III. Often authors have attended to different con- 
siderations and constraints, so the parameter k has been 
taken to be different. However, as we'll point out below, 
had the same considerations been employed, the value of 
K would have to be comparable. 

There are significant differences between the non- 
relativistic limits listed above. The most radical distinc- 
tion bundles together type I and IIB theories in opposi- 
tion to type IIA theories and the single relativistic theory 
leading to a type HI limit. Because in type I and IIB the- 
ories non-relativistic particles are sensitive to two fields, 
which mimic each other in the Newtonian regime, the 
gravitational constant is effectively renormalized. In the 
Newtonian regime (non-relativistic limit, high total New- 
tonian force), we have // « 1 or z^ sa 1, and so cj) becomes 
proportional to the Newtonian potential: 
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i.e. (f) doesn't vanish but rather shadows $Ar multiplied 
by k/Att. This "renormalizes" the gravitational constant: 
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and Giien is the gravitational constant measured, say, 
by the Cavendish experiment. Nevertheless cosmology 
(for example, Friedmann's equations) is sensitive to the 
"bare" G. Constraints arising from Big Bang nucleosyn- 
thesis therefore require k to be of the order of ~ 0.01 or 
smaller. Structure formation considerations may fix fur- 
ther the value of k (see [Sj and references therein). The 
conclusion is that in the non-relativistic regime the field 
(j) must be suppressed when qn — |V$7v| is very large. 

However, astrophysical applications of type I and IIB 
theories require that when a^ < oq the total $ must 
have MONDian behaviour. This requires simultaneously 
that 4> be in the MONDian regime and that </> be the 
dominant contribution. But this means that we must 
switch on MONDian behavior in (j) at Newtonian accel- 
erations aN much higher than oq. Only thus may the 
relative importance of 4> start increasing with decreasing 
ajsr so that when qn drops below ao it has caught up 
with $jv. Assuming the free function turns from 1 to a 
single power-law (and ignoring the MONDian magnetic 
field where appropriate) we have F^/ Fj^ ex l/\/7jv once 
MONDian behavior in (j) has been triggered. Given (|9]) 
we should therefore trigger MONDian behavior in (p for: 
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(with CjY ^10 ms for typical n) or equivalently 
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FIG. 1: Log plot of ratio between the MONDian and New- 
tonian forces, F^/Fn, against z = {k/A-n)\F^\/ao (bottom 
axis) and Fn /clq (top axis). So that Fn ~ F^ when F^ ~ ao 
(and so 2 = k/Ati; also Fn ~ ao) and at the same time 
have F^/Fn ~ n/in <C 1 in the Newtonian regime (z 3> 1, 
_Fjv — > oo), we must trigger MONDian behaviour in at ac- 
celerations much larger than ao (when 2 ~ 1). 



also much larger than ao. This point is illustrated in 

Fig.m 

This simple argument fails should fi be divergent (with 
part of its domain excised). Then, Fc/, goes to a con- 
stant as Fn — > 00, and so Gjien — G (see [l6|, Il7|). 
Consequently, it is possible to have a^^^ ~ oq, without 
fine-tuning the free-function or inducing unduly different 
GRen, in these theories. Such functions, however, may 
have other problems, rendering them non-viable. (This 
is explained further in in Sections FVl and FVI A[ see Fig.[T4l 
which contains the counterpart of Fig [1] for these theo- 
ries.) For the same reason the simple argument laid down 
in the previous two paragraphs also fails for type IIA the- 
ories, for which G is not renormalized. 

Excluding the very last proviso, our considerations 
apply equally to type I and IIB theories. We have 
parametrized the free-function v for type IIB theories 
in a way (at odds with the literature) which allows com- 
parison of "like with like" with type I theories. Thus, for 
the same k both types of theory renormalize the gravi- 
tational constant by the same amount. They also then 
predict ~ $jv for ajv ^ oq when their functions ^ 01 v 
trigger MONDian behavior at z '^ 1 and i) ~ 1, for type 
I and IIB theories, respectively. 

The features we have highlighted explain the large size 
of the bubbles around the saddle inside which type I 
and IIB theories display anomalously large tidal stresses. 
These bubbles are large (of the order rp ~ 380 km for 
the Earth-Sun saddle) because they represent the region 



where the field has started to behave in a MONDian 
fashion. This is given by the region where oat < a^^^ and 
not where ajv < ao, as might be naively expected (and 



indeed it can be easily computed that aN 



trig 
T-N 



around 



ro). It is important to stress that in a LPF saddle test we 
are probing the regime where (j) has gone fully MONDian 
but hasn't yet dominated ^n, something that happens 
at total Newtonian forces in the range ao < Qn < o,fq^^ ■ 
In spite of the dominance of $jv the MONDian signal in 
(f) can be detected because, as we shall see, it has a dis- 
tinctive spatial variation, whereas the Newtonian tidal 
stress is just a DC component. The experiment is sen- 
sitive to the time Fourier transform of the signal with a 
sensitivity that peaks at the MOND frequency (and is 
very poor for a DC component due to 1// noise). In 
contrast, in type IIA and III theories MONDian effects 
are only triggered for a^ ^ clq ^ 10^^'' ni s^^, resulting 
in very small bubbles (a few meters across). 

Even though these considerations place type I and IIB 
theories in the same basket with regards to a saddle test, 
they have a significant difference. A well known tech- 
nicality is that type I and III theories have a curl term 
(sometimes dubbed a "magnetic field"), in the sense that 
if one attempts to linearize their non- linear Poisson equa- 
tion by introducing an auxiliary vector field (e.g. /zV0 
for type I theories) this field has non-zero curl. The same 
doesn't happen for type II theories, which are already lin- 
ear in (j) (but driven by a function of the Newtonian field, 
vV^N, which has a curl). This turns out to have a sig- 
nificant quantitative (but not qualitative) effect upon the 
magnitude of the saddle tidal stresses. The magnetic field 
is known to soften the anomalous tidal stresses around 
the saddle points in type I theories, as explained in [11 1. 
As we'll show in (2a], type IIB theories have a quantita- 
tively stronger saddle signal than type I (once their k are 
adjusted to produce the same physical properties), due 
to the absence of this curl fleld^ . 

In this paper we focus on Type I MONDian theories, 
but in the conclusions explain why our results are quali- 
tatively applicable to Type IIB theories too (indeed the 
SNR forecast here are higher for type IIB theories [26|). 
Type IIA and III theories are the only ones to fall through 
the LISA Pathfinder net. 



III. THE SIGNAL TO NOISE EXPECTED 
FROM A SADDLE FLYBY 

The quantitative predictions for type I theories have 
been extensively studied using both analytical methods 
resorting to simplifying assumptions [11] and numeri- 
cal techniques [l3| including all the complications of the 



^ Note that if one defines the curl field as SJ<f> — i/V't'jv in type 
II theories (in analogy to /iV0 — V"t>jv in type I theories), then 
these theories do have a curl field. 
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FIG. 2: The transverse MOND stress signal Syy along the 
lines y — 25, 100 and 400 km (top to bottom), for the Sun- 
Earth saddle, taking the effect of the Moon into account. The 
different lines represent different lunar phases: new Moon 
(thick, black, solid) , full Moon (thick, black, dashed) and with 
the Moon appears 18° away from the Sun towards positive y 
(thin, black, solid). We also show, for the y = 25 km case, 
the Newtonian stresses (grey) rescaled by k/4-k (see text for 
discussion) . 



problem, such as the perturbing effect of the Moon and 
planets. Figure [2] has been borrowed from [Ij] to illus- 
trate the expected tidal stress along lines missing the 
saddle at 25, 100 and 400 km. As in [ij], we adopted 
a coordinate system with x aligned along the Sun-Earth 
axis and centered at the saddle and considered trajecto- 
ries parallel to x {y = b lines, where b is the impact pa- 
rameter), but other trajectories are easy to implement. 
Due to a number of practical issues |14| , only transverse 
tidal stresses can be measured, say the Syy component. 

Predictions are cast in the form of tidal stresses be- 
cause this is what is directly measured by the instrument. 
LPF measures the relative distance between the masses, 
but this is converted into a relative acceleration (or its 
Fourier transform in time) . Up to a factor dependent on 
the proof mass separation, the measurement is therefore 
one of tidal stress along the direction linking the two 
masses (with further masses, other tidal stress compo- 
nents would become accessible). In line with this state- 
ment, noise evaluations and forecasts are expressed in 
terms of tidal stress or relative accelerations; one should 
use the inter-mass separation to convert between the two. 

It is of paramount importance to note that the field (j) 
produces both a MONDian effect and a Newtonian pat- 
tern, associated with a rescaling of G in the Newtonian 
limit. The properly MONDian stress is therefore: 
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dxidxj At: dxidxj 



(13) 



i.e. we must subtract from (f> its component included in 
the Newtonian background, which is $Ar rescaled by n/in 
(see [i3| for more details). In performing this exercise 
it is essential that and ^^ have been solved to the 
same degree of accuracy. In Section IIVCI we will discuss 
the impact of an imperfect subtraction of the Newtonian 
component. 

The data analysis task in hand is therefore to detect 
a "wave form" of this type with the instrument aboard 
LPF. As a first hack at the problem, we evaluate the 
performance of noise matched filters. Matched filter- 
ing is a well-known data analysis technique used for effi- 
cientl y d i ggin g a signal with a known shape out of noisy 
data [1^ |23|. The technique is extensively used in the 
search for gravitational waves. The idea is to correlate a 
time series x{t) with an optimized template designed to 
provide maximal signal to noise ratio (SNR), given the 
signal shape h{t) and the noise properties of the instru- 
ment. Generally we have x(t) = h{t — ta)+n(t), where ta 
is the signal "arrival time" and n(t) is a noise realization. 
We want to correlate x{t) and an optimal template q(t), 
yet to be defined, according to: 



c{t) 



x{t)q{t + T)dt 



(14) 



where r ia a lag parameter, giving us essential leverage 
if we don't know ta a priori. The average of c over noise 
realizations is the expected signal, S, and its variance is 
the square of the noise in the correlator, A^^; the forecast 
signal to noise ratio is therefore p = S/N . A straight- 
forward calculation (under general assumptions, namely 
the Gaussianity of the noise — more on this later) shows 
that p is maximized by choosing a template with Fourier 
transform: 



q{f) 



qit)e'^^f^dt = ^^^^; ... (15) 



and setting the lag t to the arrival time, t = ta- Here 
Sh{f) is the power spectral density (PSD) of the noise, 
conventionally defined from 



{h{f)n*{f)) = ::S,{f)S{f~r), 



(16) 



(the factor of 1/2 hails from the tradition of taking one- 
sided Fourier transforms of the noise auto-correlation — 
-i.e. with / > only). The maximal SNR, realized by 
the optimal template, is then: 



Popt 
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(17) 



Notice that the optimal template, q{t), defined by (fTS)) 
is not the signal, h{t), but rather a filtered version with 
a pass where the noise is low and a cut where the noise 
in high. Also, the optimal SNR given by pT|) is not 



the energy in the signal but an integrated signal power 
weighted down by the noise PSD. 

These techniques are run of the mill in gravitational 
wave detection, where the arrival time of a signal is often 
not known^. For example, for a chirping signal, even if we 
have a fair idea of the shape of the signal, we can't know 
when a binary coalescence is to take place. We therefore 
have to shift the template Fourier transforms, h{f), by 
all possible phases, until the maximal SNR is obtained, 
should there be indeed a signal. This adds an extra pa- 
rameter to the fit and may also be the source of spurious 
detections. It affects the management of 1// noise and 
increases the false alarm rates (as effectively we have a 
number of trials equal to the total observation time di- 
vided by the duration of the template). This problem is 
absent in the context of our test, where ta is known since 
we do know where the saddle is and therefore where the 
signal is meant to start in the time-ordered series. A 
natural truncation in integration time T is also present, 
simplifying 1// dealings. 

It has been estimated'^ that the saddle can be pin 
pointed to about a kilometer and the spacecraft loca- 
tion determined to about 10 km even with most basic 
tracking methods*. The effect this has on SNRs is neg- 
ligible, indeed the SNR grids we are about to present 
have this sort of resolution. Thus, we can simply set 
ta = with an appropriate choice of conventions and set 
to zero the time lag t in the correlator c, to achieve opti- 
mal results. For all practical purposes the starting time 
is indeed known. To the same degree of approximation, 
we also know the spacecraft trajectory and velocity with 
respect to the saddle. 

Given a spacecraft trajectory, the conversion of tidal 
stresses (such as those depicted in Fig. [2]) into a template 
in time, h{t), is then trivial. For a setup such as the one 
described above we have h{t) = Syy{vt,b,0), where v is 
the velocity of the spacecraft, and t — corresponds to 
the point of closest saddle approach. In a more general 
setup, for an approximately constant velocity v, a closest 
approach vector b, and with the masses aligned along 
unit vector n, we have 



h{t) =n'n^S,j{h + vt) 



(18) 




FIG. 3: The amplitude spectral density (ASD) of the MOND 
tidal stress signal for a trajectory with 6 = 50 km and 
V = 1.5 km s~^, compared to the ASD of the basic noise 
model described in the text, assuming a baseline of 1.5 x 
10"^" s"^/VHi. This scenario generates a SNR of 28. 




This template should be Fourier transformed and, given a 
noise model, used to produce an optimal template, using 
a noise matched filter. Its SNR can then be evaluated. 
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FIG. 4: Signal to Noise ratio contours, for various impact 
parameters up to 600km and base noise ASD. We set the 
spacecraft velocity at 1.5km s~^. Calamitous assumptions 
would still lead to SNR of 5. More optimistic ones (6 around 
50km, noise half way up the scale) would lead to SNRs easily 
around 50. 



There are exceptions, for example if the signal comes from a 
supernova or any other source for which there is an extrinsic 
method, typically in the optical domain, for flagging the source 
of gravity waves. 

These uncertainties are of a practical nature and should not be 
confused with theoretical uncertainties. It can be estimated that 
the MOND saddle cannot be shifted with respect to the New- 
tonian saddle by more than a meter, and this is just an upper 
bound. 
S. Kemble, private communication; to be published. 



To gain some intuition on the nature of the signal in 
Figure |3] we plot the amplitude spectral density (ASD) 
of the signal, which is the square root of the PSD: 



P{f) = ^ 



+T/2 



T/2 



dt h{t) e-^""'^' 



(19) 



where / is the frequency, t is the time and T is the in- 
tegration period (here taken conservatively to be T = 
2 X 10* s). This can be directly compared to the noise 



ASD, the form usually quoted by experimentalists. As 
a simplified LPF noise model (see [ij]), we assume that 
the noise is white in the frequency range between f and 
10 mHz, i.e. we assume a constant baseline with ASD 
around 1.5x10"^"' s~^/-\/Hz. For lower frequencies we as- 
sume 1// noise and for higher frequencies that the noise 
degrades as f^. With these assumptions the noise and 
signal ASDs are plotted in Figure |31 for typical param- 
eters. As we can see, there's signal to noise of order 10 
over a couple of decades, making it not surprising that the 
integrated SNR is in double figures (in this case around 
28). 

We can now run through the parameter space of the ex- 
periment and evaluate SNRs. For example, let's assume 
ti = 1.5 km s~^ and explore impact parameters up to 600 
km. Let's also consider the effect of changing the base 
line ASD of our noise model. The result is plotted in Fig- 
ured We see that we'd need to miss the saddle by more 
than 300 km to enter single figures in SNR, with typical 
noise levels. For b around 50 km a SNR of the order 30-40 
is not unrealistic. Recent work has placed a figure on the 
impact parameter around 6 = 10 — 50 km within reach. 
In combination with the expectations for the noise, this 
makes the test very promising indeed. However we should 
now look at this preliminary analysis in more detail. 




FIG. 5: Plot of SNR against satellite velocity for an im- 
pact parameter of 50 km and a baseline noise of 1.5 x 
s~'^/vHz. We note a broad peak around v = 2 km s~^. 



10" 

Higher speeds shift the signal to higher temporal frequencies; 

however the rough speeds of all trajectories in the Earth-Moon 

system are already optimal, given the noise properties of the 

instrument. 



IV. FURTHER DISCUSSION 

In this Section we refine and discuss further the basic 
results presented in the previous Section. There is con- 
siderable uncertainty regarding the details of the flyby 
trajectory, namely its speed. In Section IIV Al we show 
that the effect of the speed is minimal, within the range of 
speeds expected from any trajectory in the Moon-Earth 
system. In Section IIVBI we present improved, more re- 
alistic noise models, repeating the analysis with a best 
and worst case scenario for instrument performance as 
understood at the time of writing. We also outline work 
in progress, improving on noise matched filters and on 
estimates of false alarm rates. Finally in Section llVCI we 
discuss issues related to the background tidal stresses, 
namely the the Newtonian background and the space- 
craft self-gravity. 



A. The impact of the spacecraft velocity 

What is the effect of the spacecraft velocity on the 
SNRs presented in the previous Section? The question 
is relevant as it can assist the strategy in designing flyby 
trajectories. Here we show that in practice all that mat- 
ters is the trajectory location (impact parameter and pos- 
sibly angle). Within the range of realistic speeds, the 
SNRs do not vary substantially. The good news is that 
due to a remarkable coincidence, these speeds are already 
near optimal. 



As Eq. (fT5)) shows, the spacecraft velocity is the con- 
version factor between the spatial scale of the tidal stress 
and the time scale at which the instrument measures 
them. Of course, in detail, this has an effect on expected 
SNRs. Higher/lower speeds mean a faster/slower scan- 
ning oj these spatial features, and thus a shift of the tem- 
plate h{f) to higher/lower frequencies, whilst keeping the 
noise ASD fixed. Therefore the SNR has to change. This 
is shown in Fig. [S] for 6 = 50 km and a baseline noise of 
1.5 X 10-14 s-VVH^. 

We see that the SNR has a peak dX v — 2 km s~^ . 
However this peak is very broad with respect to the type 
of variations that might be expected from different tra- 
jectories leading from LI to the saddle [Tj]. For the rough 
range v = 1.5 — 2.5 km s"^ the SNR varies in the range 
27-28, approximately. For w = 1 — 3 km s""'^ (which is 
pushing it, in terms of real orbits) the variations would 
be in the approximate range 25-28. The message is clear: 
get as close to the saddle as possible, never mind the 
speed. The speed will never be very far off the optimal. 

This result can be understood qualitatively. As a crude 
estimate, anything moving in the Earth-Moon system has 
a typical speed of the order of 1 km s^^. The MONDian 
tidal stress for the Earth-Sun saddle displays variations 
on a length scale of the order of 100 km. Therefore the 
MONDian signal will always be felt by LPF on a time 
scale of minutes, i.e. in the mHz range. This is just 
where the instrument noise is lowest, a remarkable co- 
incidence considering that the instrument was built to 
these specifications for entirely different reasons (astro- 
physically motivated gravitational wave templates have 
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FIG. 6: This figure replots Fig. [S] adding on the best and 
worst case scenarios for more realistic noise models, as at 
the time of writing. We have assumed a trajectory with the 
geometry described in the main text, with impact parameter 
of b = 50 km and velocity v — 1.5 km s~^. We have also 
plotted the contribution of to the Newtonian background. 



these time scales) . And yet the typical speeds and length 
scales of the problem combine to make the instrument al- 
ready optimal for a MONDian saddle test. 



B. Improved noise models 

A number of improvements to the noise model used in 
the previous Section are possible. These are the subject 
of a paper in preparation [28|. Obviously there isn't a 
frequency region with white noise. Instead, the noise is 
likely to be higher than modeled in Section IIIII at high 
frequencies but lower than expected at low frequencies. 
The turnover between the two regimes is smooth, as de- 
picted in Figure ini where we superimposed the simplified 
noise model used in Section IIIII with the more realistic 
estimates for ASD for a best and worst case scenario. It 
has been argued that the worst case scenario might be 
too pessimistic and the best case scenario too optimistic, 
so we should take these two models as extremes. 

In Figure[7]we plot the SNR as a function of impact pa- 
rameter with u = 1.5 km s~^, assuming the two extreme 
scenarios. As we can see, in the best case scenario we'd 
need to miss the saddle by more than 650 km for the SNR 
to drop below 5. In the worst case noise scenario, how- 
ever, that figure would shrink to about 250 km. For the 
currently expected & ^ 50 km the SNR would be in the 
range 13 — 44. In spite of the uncertainties, all scenarios 
lead to optimistic prospects (and even overkill) regarding 
a detection. We stress that we will know what the noise 
is, in situ and while on LI. Our forecasts are useful, but 
we should highlight that they'd become concrete, fixed 



FIG. 7: The SNR for the improved noise models (best and 
worst case scenario) assuming t; = 1.5 km s~^ for a variety of 
impact parameters b. 



numbers once the mission goes ahead. 

We should add that even if the noise ASD is known, 
further issues complicate the simple data analysis proce- 
dures presented in the previous Section. Most notably 
the real noise is non-Gaussian and non-stationary. This 
may increase the probability of a "false alarm" , to use the 
jargon of gravitational wave detection. Putting a realistic 
figure to the probability of a false detection requires hav- 
ing the instrument switched on before and after a saddle 
flyby, characterizing the noise in situ, and evaluating the 
false alarm rates with real noise. No prior modeling can 
be a substitute for this. Nonetheless more realistic simu- 
lations of the instrument response and noise are possible. 
We are currently working on these. 

The issue of false alarm rates is obviously central 
should there be detection. But even just planning the ex- 
periment, it raises important questions, e.g.: given these 
rates, is it better to sacrifice b at the expense of multi- 
ple flybys, or should we put all our efforts into a single 
flyby with a 6 as low as possible? Should the noise be 
approximately Gaussian and stationary, the probability 
of a false detection is simply [l^ : 



T = iVcrfc(p) 



(20) 



where p is the optimal SNR, and N is the number of 
trials. In gravitational wave detection N = O/T, where 
O is the total observation time and T the useful dura- 
tion of the filtered template. The factor N can be very 
large, so that even substantial SNRs (say 8 or 9) can 
produce non-negligible rates J^. In gravitational wave 
detection this nuisance can be mitigated by coincident 
observations. We stress that no such problem is present 
here. We do know where the saddle is for all practical 
purposes, so N = 1, removing the extra factor enhancing 
the false alarm rate. 
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FIG. 8: ASD plot of the MONDian and Newtonian signal 
(multiplied by n/An), as compared to the noise ASD. We con- 
sider the effect of subtracting the Newtonian component in 0. 
This only affects very small and very large frequencies. 



The high SNRs we've obtained at low b suggest that 
it would not be advisable to sacrifice b for the sake of 
multiple fiybys, in order to reduce false alarm rate. This 
statement should be further scrutinized using real noise. 
But even if it's true there is an important sociological el- 
ement. The reliability of any scientific claim rests on re- 
producibility. Should there be a positive detection, more 
than one flyby would go some way towards establishing 
the case for reproducibility. 




C. The New^tonian background and self-gravity 

We finish this section by examining two possible sys- 
tematics that could plague a saddle test: the Newtonian 
background and the spacecraft self-gravity. These are 
natural concerns, but their impact is negligible. In estab- 
lishing this fact it is important not to confuse force and 
tidal stress. It is also essential to examine the Fourier 
components of the stress signal and distinguish a DC 
component from a signal peaking at frequencies to which 
the experiment is sensitive. 

The MONDian saddle signal has a spatial scale tq « 
380 km. In this region, apart from an inner bubble a 
few meters across, the Newtonian force is always much 
larger than the MOND force and also the MONDian ac- 
celeration ao ~ 10^^" m s^^. We recall the discussion 
in Section ini where we noted that LPF would probe the 



< a 



trig 



N 



with 



trig 



'AT 



In- 



regime ao < un 
deed around tq ~ 400 km, the Newtonian acceleration 
is ttM ^ o,N^ ■ "^^^ Newtonian tidal stress is therefore 
dominant in this regime (with an intensity of the order 
A ^ 10~^^ s~^), bu t, crucially, it is approximately a DC 
component [ll|, [ij] . This is to be contrasted with the dis- 



FIG. 9: This plot illustrates the systematic effects that might 
result from an incorrect Newtonian subtraction. We consider 
the transverse tidal stresses felt in trajectories with impact 
parameters 6 — f 00, 500, 1000 km. We then subtract the DC 
constant Newtonian tidal stress contributions from (j) (top) 
and its full contribution (bottom). As we can see an imperfect 
subtraction produces a spurious ramp in the stress. 



tinctively varying MONDian signal (see Fig. [2]) which, as 
we've shown, translates into a signal peaking at frequen- 
cies where the noise is low. A DC component, on the 
other hand, is well buried in the 1// noise. It is true 
that in detail the Newtonian tidal stress is not exactly 
constant on the scale of tq. But we do know what it is, 
to the same accuracy as we know the saddle location and 
trajectory, and can subtract it off. Furthermore its spec- 
tral shape away from its DC component is very different 
from that predicted by the MOND signal, as shown in 
Fig. m The Newtonian background amounts to the sub- 
traction of a known component. 

A related matter was flagged in Section IIIIl and re- 
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lates to subtracting off from its effect on the renor- 
malization of the gravitational constant. Some of con- 
tributes to the Newtonian background and should not be 
included in the MONDian predictions (cf. Eqn. P^ : see 
also [13]). In Figure [9] we plotted the power spectrum 
of the non-DC component of the Newtonian tidal stress 
produced by 0. The impact of not subtracting the com- 
ponent of (f) contributing to the Newtonian measurement 
can be appreciated in Figure ID This figure also gives us 
an idea of the level of impact an imperfect Newtonian 
subtraction might have. We considered the transverse 
tidal stresses felt in trajectories with impact parameters 
b = 100, 500, 1000 km. In the top we subtracted only the 
DC component, in the bottom we subtracted the full con- 
tribution of (j) to the Newtonian tidal stresses. As we can 
see an imperfect subtraction would produce a spurious 
ramp in the stresses. 

Another issue is LPF's self-gravity. The mission re- 
quirement is that the differential acceleration of the two 
LPF test masses should be balanced at the level of 
a ~ 10"^ m s~^, but actual performance may beat the 
nominal requirement by a factor of 10. Yet again this is 
a DC component and does not affect the measurement 
in tidal stresses with the distinctive temporal variations 
we have posited. There are of course time- varying uncer- 
tainties in the self-gravity balancing but these are much 
smaller. They are mainly due to thermoelastic effects. 
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Hz at least 



and are on the level of 3 x 10" 
down to ImHz. 

An issue related to this concerns the position of the 
saddle. Naively one might think that with a self-gravity 
of the order of 10~^ m s~^ the position of the saddle 
would be perturbed by the spacecraft. The two test 
masses could even generate distinct saddle points due 
to their gravity. This concern ignores the fact that with 
realistic impact parameters we are not testing the regime 



UN 



ao, but the regime oq < un < 



trig 



with much 



larger Newtonian accelerations. For instance for an im- 
pact parameter of 40 Km we have ajv ^ 10" 



Around b 



trig 



ro 



the Newtonian acceleration is 



ms 



10- 



We'd need to approach the sad- 



dle much closer than about 400 meters before self-gravity 
becomes an issue and the spacecraft itself had to be in- 
cluded in the computation of the location of the saddle. 



V. MONDIAN FREE-FUNCTIONS 

In the second part of this paper we examine the gen- 
erality of our predictions. So far we have focused on 
type I theories with a specific fitting function /i (the one 
used in |13|). But even if we restrict ourselves to type I 
theories there is a whole free function ii{z) to play with. 
Would theorists be able to wriggle out of a negative result 
availing themselves of this freedom? In the next 3 sec- 
tions we prove that under general conditions only type I 
theories with fine-tuned /i-functions would survive a neg- 
ative result (with notable exceptions to the rule, spelled 



out). This conclusion is also expected to apply to type 
IIB theories, although we won't prove it in detail. Type 
III and type IIA theories turn out to be the only ones to 
evade a LPF saddle test. 

In this Section we start by reviewing previously pro- 
posed fj,, laying down a common notation. We then dis- 
cuss criteria for physically permissible /i, showing that 
for single power-laws they all produce SNRs of the same 
order, for impact parameters smaller than 400 km. In 
Section I VI I we will then consider more structured func- 
tions fj,, with more regimes and scales. 



A. Notation and previous proposals 

As explained in Section [TTl for type I theories two po- 
tentials act on non-relativistic test masses: the Newto- 
nian potential ^^v and a fifth force 4>. Thus, the total po- 
tential is $ = $Ar -I- 4>, or in terms of forces, F = F^r -|- F^ 
(as before, we set S = 1 for simplicity.). We recall that 
both contributions satisfy Poisson type equations: 



V^$ 



N 



AttGp 



V • (AiV0) = nGp 



(21) 



(22) 



where, as before, we write u with argument z = ■^- — —. 
In the Newtonian limit /i — ?> 1, whereas MONDian be- 
haviour in (f) is triggered when fi ^ z. Note that here 



z = 



(23) 



where y is the variable eniployed by Bekenstein in his 
original paper on TeVeS [7]. Much confusion has arisen 
from different notations in this respect. 

We should not confuse /i(z) with the function fl{x) 
used in type III theories (cf. Eqn. ([7])) and also favoured 
by astronomers. Even in type I (and also II) theories, we 
can loosely define an effective fl{x), obtained from adding 
equations ([2T|) and (|22t and comparing with Eq. ([T]). 
This effective fl{x) function is frequently used in fits to 
galactic phenomenology. However the two functions jl{x) 
and /i(z) can only be easily related if the MONDian curl 
term can be neglected. This proviso is often incorrectly 
ignored. If the curl term is non- negligible, then type I 
theories don't properly have a il{x) function, and there's 
no substitute for integrating the equations on a case by 
case basis. 

If the curl field can indeed be ignored in the integration 
of (p2|) . then it's easy to relate functions fj,{z) and fl{x) 
(see, e.g. Q). Using ^Ul and (f^ . their definitions can 
then be rewritten as F = Fn/JI and F0 = -^Fj^, so 



that F 



■ N 



Fa, implies: 



(24) 
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In addition we can write the argument x — F/ag in terms 
oi z = {K/ATT)F^/ao by deriving: 



47r 
X — — z I 1 

K 



4:Trfi{z) 



(25) 



Eqns. ([M]) and ([25]) provide a parametric expression for 
/i(x). Note that Eq. (IMl) triviahy imphes that in the 
Newtonian regime (/i « 1) MOND has the effect of renor- 
mahzing G as: 



GRen = — ~ G{l + -—] 



(26) 



a result aheady presented in Section |TT1 

Several /i functions have been previously proposed. 
The "toy" model used in Bekenstein's original paper Q 
follows from the implicit expression: 



,_^lHu-2r 



(27) 



4(1 -m) 

A variation was employed in [ll| to facilitate analytical 
work on the 2-body problem: 



z2 = 



/r 



(28) 



1-^4 

In some manipulations an inversion of the latter is useful 



M 



Vl + 4z4 - 1 

2^2 



(29) 



A proposal quite distinct from these two can be found 



^(z) 



1 - 



(30) 



with the case a — I first suggested in [2l| . To bridge our 
notation with the fJ.s{s) used in [l^ we should use the 
dictionary (obtained from direct comparison of (j2ip and 
(|22|) and their counterparts in [iq): 



M 

z 



i^'^^ 



47r 



(31) 
(32) 



We stress that this function diverges, a property that sets 
it apart from all those functions which tend to a constant 
as z — > oo. Underlying this statement is the postulate 
that the domain of the function should be excised after 
the divergence is reached, i.e. one should impose s < 
1/a. The distinction between bounded and divergent 
/i was extensively studied in [17] and is indeed central 
to the discussion^. A hybrid possibility, incorporating 



the behaviour of (|30p on galactic scales into a bounded 
function, can be adapted from the proposal in |2l| . and 
will be examined in Section IVI Al 

There is some debate over which /2 functions best fit 
astrophysical data. Examples include [29i] : 



and 

or even [H, [U 

/2(.t) = 



fl{x) 



jj.{x) 



VT 



2x 



1 + (2 - a)x + yj{l - axf + 4a; 



(33) 



(34) 



(35) 



With the proviso spelled out above (non-invertibility of 
a /i(z) in terms of a /i(x) in the presence of a curl field), 
function ([55)) can be derived from ([50]) . Likewise (P7)) 
and (HSl) lead to: 



li.{x) 



1- Vl + 4a: 
2x 



^/YT^i-\ 



Vl + 4a; + l 
2x 



1 + 2a; + vT+4x 



(36) 



where we have written three algebraically equivalent ex- 
pressions to facilitate comparison with the literature. 
Note that although (|36p follows from ([35]) for a = 0, the 
same doesn't happen with their /i functions, and Beken- 
stein's proposal (P7)) is strictly not covered by ([50]). The 
claim has been made [13, llll that galactic observations 
favour a — \^ . 



B. Permissible, non-fine tuned [i functions 

Putting aside detailed predictions for galaxy rotation 
curves (which may well have been combined with incon- 
sistent approximations, e.g. regarding the curl field), the 
following criteria are reasonable for physically permissi- 
ble, non-fine tuned /i functions defining type I theories: 

• A. The cosmologically measured G cannot dif- 
fer significantly from that measured, say, by the 
Cavendish experiment. That is: Grew ~ G. 



^ As a notational word of caution note that for convergent free- 
functions we have defined \i, such that ^ — > 1 as 2 — > oo; whereas 



in [l7j one has /^s 

^0 = 47r/K. 



/iOj so that in effect we get the dictionary 



^ As a further notational word of caution (c.f. llTll ) note that our 
jl tends to G/Gn^n (or in th^ notation of [ITII . to I/vq). The 
Milgrom-like proposals considered above all tend to 1. Thus, 
these proposals can only be approximately true with Guen ~ G. 
Otherwise, the proposals considered here should strictly speaking 



be labelled fi.^''^'-3{x), with /i : 
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• B. When the total Newtonian acceleration ajv 
drops below uq the full potential $ must be in the 
MONDian regime, that is, we need (p to be in the 
MONDian regime and to dominate <i>jv. 

• C. Function ^ should only have one scale, below 
which (j) is MONDian, and above which it is near 
Newtonian. The detailed form of the transition is 
left undefined, but fi should have a single transition 
from 1 to z. 

Items A and B have already been discussed in Section HIl 
Item B is the most basic requirement for the theory to 
be of astrophysical use, regardless of the details. Item C 
has been spelled out because it will be broken in the next 
Section, to illustrate just how finely tuned /i would have 
to be to evade a negative saddle result. 

As explained in Section |lTl these requirements imply 
that (j) must enter the MONDian regime at a much higher 
acceleration than oq, leading to an intermediate regime 
Go < ttN < a^^ where (f> is fully MONDian but still sub- 
dominant to $Af. This implies that for any fi satisfying 
these constraints, when un '^ ao (i-C- for astrophysical 
applications) we must necessarily have 




200 300 400 

b/km 



600 



FIG. 10: Signal to Noise ratio contours, for various impact 
parameters up to 600km and base noise ASD, using the same 
templates as in Fig. 3] but with an exponential fall off in the 
model-dependent region r > ro. As we can see, for impact 
parameters b > 400 km the SNR drops more sharply, but 
nothing changes very much for b < 400 km. 



F^ w V FArao 



(37) 



This statement is independent of k and only relies on the 
fact that n M z ^ :^^-^ in the MONDian regime. // 
the curl term can be ignored we therefore have zF^ = 
-^Fm, and thus (|37p follows. Recalling x = F/a^ we 
must conclude that: 



III theories). The predictions for hit) for r < rg are also 
model independent, since they rely on /i w z, for z < 1. 
However the predictions referring to regions with r > r^ 
depend on the exact form of the transient from /x sa z to 
/x « 1, because they depend on 5^1, not fi. For example, 
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(38) 
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(40) 



The exception to this rule is obtained with a diver- 
gent /i, as already announced in Section |IT1 Then, we 
discover the interesting behaviour |la, [T^l that F^p goes 
to a constant as Fn grows to infinity, instead of becom- 
ing proportional to Fn (c.f. Eq.®). Specifically, taking 
model ((50)) . we find that: 



^0 ~ — 
a 



(39) 



so that asymptotically no renormalization of G takes 
place: G_Ren = G. Such a functions would lead to a 
different p,, as we've seen. In Section IVI Al we'll show 
that the MONDian behaviour driven by these functions 
would be rendered invisible to LPF. 



C. SNRs and /j, dependence 



If we take the whole class of /i satisfying requirements 
I, II and III we conclude that they have the same a^^^ 
and consequently the same tq . Fiddling with fi therefore 
doesn't change the spatial scale of the effect for type I 
theories (and also for type IIB, but not type IIA or type 



has a very different fall off from 



VT 



1- 



2z2 



(41) 



a point recognized in (ll[. Missing the saddle by a lot 



more than 400 km would therefore leave us at the mercy 
of model dependence, and fi functions satisfying A, B, C 
could be found bypassing a negative result, e.g. 



M = 



(l + zf^)^"' Pz^ ^^'^ 



with a large /3. However, for trajectories hitting the re- 
gion r < To (i.e. for h < 380 km) the peak of the signal 
is actually model independent, and therefore the SNRs 
predicted aren't expected to depend on the details of the 
theory. 

As an extreme illustration of the model (in) dependence 
of our SNR predictions we have excised the signal out- 
side the MOND bubble from our templates, imposing an 
exponential fall off. Fig. [TU] is the resulting counterpart 
to Fig. 01 We see that for b < 400 km our conclusions 
remain substantially the same. For b > 400 km the SNRs 
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drop much more sharply. This is the worst that can be 
expected. 

Impact parameters around 50 km are now considered 
within easy reach. In order to bypass a negative result we 
would therefore have to shrink the bubble size, defined 
by scale tq. This requires breaking condition C and con- 
sider contrived fi functions with two scales, which we now 
proceed to do in order to appreciate the full implication 
of a negative result. 



VI. A NULL RESULT AND DESIGNER ^ 
FUNCTIONS 



It is often difficult to falsify a theory containing free pa- 
rameters: all that can be readily done is to constrain its 
parameters. However the constraints may be such that 
the theory becomes contrived beyond "reasonable". In 
what follows we imagine a scenario where no anomalies 
are found with respect to the Newtonian expectation, for 
b < 400 km. Obviously all the theories considered so far 
would be ruled out, to a degree of significance of the same 
order as their expected SNR. The issue would then be- 
come to determine which "designer" functions /z predict a 
SNR of order 1, thereby surviving a "no anomaly" result. 
The more contrived the required fi, the more blatantly 
one should throw in the towel. 

In proposing a designer /i we shall impose that it sat- 
isfies requirements A and B to the same extent as the 
functions we've been considering. The theory should 
still be of astrophysical use and not confiict with ob- 
servations on very general grounds. However, we drop 
requirement C, allowing the function to have two inde- 



pendent scales (notice that 



trig 



is not independent for 



the models considered so far). Specifically, we endow fi 
with an intermediate power n ^ 1 linking the Newtonian 
regime, where /x « 1, and the astrophysically relevant 
MONDian regime, where fx~z. 

Requirement B demands that jjl ~ z for z < k/Att, 
as before, so that F^ w Fn when Fn w aq, and F^ w 
^/FnUq for ttN < oq- Requirement A imposes /i — > 1 
for large z, so that Gren is the same as for the single 
power-law /i considered before (cf. Eqn. ([25)) ). If we are 
to shrink the size of the MOND bubble so as to accom- 
modate a negative outcome from a saddle test, then we 
need a sharper power, n > 1, bridging these two regimes. 
Thus F^/Fm could increase faster, with decreasing an, 
from its small value K/47r in the Newtonian regime, to 1 
at a^ — flQ. This would reduce a^^^ and thus rg. The 
point is illustrated in Fig. (|lip . where we have replotted 
Fig. [T] (made for a /i with a single power-law; recall the 
argument in Section [II| . 

These considerations fully specify the function /i, up 
to details on the transition regions. We should have: 



fL 



\ztrigj 



for 
for 



z < 



in 



47r 



< z < z 



trig 



(43) 
(44) 
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FIG. 11: Log plot of ratio between the MONDian and New- 
tonian forces, F^/Fn, against z = (fc/47r)|F^|/ao (bottom 
axis) and Fn /aa (top axis). So that Fn ~ F^ when F^ ~ ao 
(and so z = k/Ati; also Fn ~ tio) and at the same time 
have F^/Fm ~ k/Itt <C 1 in the Newtonian regime {z ^ 1, 
Fn — !> oo), we must trigger MONDian behaviour in (f> at 
accelerations much larger than ao. However, by allowing a 
sharper intermediate power-law in /i, the trigger acceleration 



trig 



may be smaller (in this illustration by a factor of 10). 



M 



1 



for 



z > z 



trig 



(45) 



where the point where non-Newtonian behaviour in (j) is 
triggered can be interchangeably pinpointed by: 



ytrig 



trig 



trig 



UttJ 



"">4^ 



flQ 



V47r/ 



(46) 

(47) 
(48) 



Notice that 



trig 
I'M 



is a now a truly independent parameter 



of the theory (which can be traded for n) . We still have 
that when oat < ao the field (j) dominates $Ar as per 
requirement B, but now the intermediate region, where 
(j) hasn't yet dominated but is already non-Newtonian, is 
in a narrower band of accelerations oq < a^ < o.^^ ■ As 
a result the MOND bubble shrinks according to 



ro 



Ht, 



km 



(49) 



In Figure [T2] we have plotted this dependence. As can 
be seen, it's easy to change rg by an order of magnitude 
with n not much different from 2. To reduce tq by more 
than that, however, a very extreme intermediate power 
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FIG. 12: The size of the MOND bubble as a function of the 
intermediate power n. It is easy to collapse to bubble by 
an order of magnitude (say to around 20 km) with n ~ 2. 
However, to make the bubble much smaller (say, on the order 
of a few kilometers), very dramatic intermediate powers would 
be required. 



x10 




b/km 

FIG. 13: Contours of the power n needed to obtain SNR= 1, 
for different noise levels and impact parameters up to 6 = 
600 km. For n 7^ 1 the function is "unnatural". We see that 
as soon as we plunge deep into the MOND bubble, a rather 
unnatural designer /i becomes necessary to accommodate a 
negative result. 



would be required^ 



Regrettably we can never make a model independent 
statement on what n is needed for a SNR of order 1. 
If nothing is observed then by the nature of the problem 
we must be making observations in the regime b ^ rQ{n). 
Therefore we are necessarily probing the transient from 
/i ex z" to /i ^ 1, dependent on the exact form of the 
function fi. Nonetheless it is interesting to perform this 
exercise, assuming a specific function, say: 



/i(z) 






For z ^ z*"3 this can be expanded as 
fj, ^ 1 + 6fi — 1 



(50) 



(51) 



Also for b ^ ro(?T.) the curl field is negligible, so we can 
write: 



^F0 = ^Fn 



(52) 



and solve it perturbatively. Expanding as in F^ = "^F^ 



6Frh, we have to zero order ^F^, = -vzFn- To first order 



"^ Notice that with this particular model the MOND bubble can 



never shrink smaller than -^380 km. 



we then obtain: 

K ,„ , /47r On \ , , 

.F,.--(W..(-^j F. (53) 

from which the tidal stresses can be inferred. 

The results are condensed in Fig. [131 depicting the 
value of n needed for a given b and noise level in order 
for a SNR of one to be obtained (and so a negative result 
be acceptable) . As we can see as soon as we plunge deep 
into the MOND bubble, a rather unnatural designer /i 
becomes necessary to accommodate a negative result. 



A. Motivated functions with features 

In this paper we have attempted not to mix galaxy ro- 
tation fits with our considerations. The reason is that it's 
not clear how these fits would stand if performed together 
with the need to fit Solar system data and a saddle test: 
the performance of goodness of fit statistics under joint 
constraints might be very different. Also penalization 
for extra parameters (such as a) has probably not been 
properly enforced (see [30] for an example of the implica- 
tions), and it would certainly behave very differently in a 
joint fit, where the number of degrees of freedom would 
be much larger (c.f. the behaviour of the Bayesian infor- 
mation criterion as seen in [30|). Having said this, it may 
well be that what we labeled "two-scale" or "designer" 
models is precisely what is needed for one such joint fit. 
We therefore examine how these motivated functions fare 
in terms of SNR for a saddle test. 

One type of function which evades a saddle test are 
those fj, which diverge (eg. [ig]). As explained before. 
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FIG. 14: Log-log plot of ratio between the MONDian and 
Newtonian forces, F^/Fn, against z = {k/4n)\F,i,\/ao for 
functions (|55|l . with n = 0,1,2. For n = we realize the 
divergent function (|30p and as we can see there are only 
two regimes, corresponding to two power-laws, with the ratio 
never flattening to a constant. For all other cases we have a 3 
piece function, with a fall off to the Newtonian regime which 
depends crucially on n. 



then the asymptotic G is not renormalized, since F^ tends 
to a constant, 






a Fn 



(54) 



as F/v —> oo. Thus the profile of F^/F^ is actually not 
"designer" in the sense we have been using the word, 
since it merely changes from one power-law (I/^/Fm) 
to another (1/i^Ar), never leveling off into a constant. 
The counterpart of Fig. [T] for this model is plotted in 
Fig. [TH curve labelled n ~ Q. Since F^/F^ never levels 



off, GRen = G, with a 



trig 



cq. a detailed analysis of the 



saddle effects of this theory cannot be carried out with 
the methods proposed in this Section. Since /.t doesn't go 
to a constant, strictly speaking full MOND-like effects 
are present for all accelerations, and the curl field can 
never be neglected. Still, one may expect the order of 
magnitude of the predicted effects to be small, and the 
associated saddle bubble to be invisible for LPF. We de- 
fer to a future publication a more detailed examination 
of the predictions associated with these functions. 

Divergent /i functions, however, may fall foul of So- 
lar system constraints (see discussions in |17H20| , but we 
stress that this conclusion might not be fully general). A 
compromise can be struck by combining the functional 
form of the proposed unbounded ^ with a curve flatten- 
ing to a constant beyond the scales probed by galaxy 
rotation curves. With such constructions we are plainly 
entering the realm of the "designer" (or multiple func- 
tional form) functions, as defined in this Section. Such 



multiple regime, multi-scale, functions are implemented 
in one way or another, in all such proposals |2l|-|24l| . For 
example, translating [2l| into our variables leads to: 



Attz 



A* 



1 



^ + a^l{l-^iY 



(55) 



On galactic scales this reduces to ((30)) . yet for the purpose 
of the Solar system we have: 



M«l 



aF,, 



(56) 



The F^/Fm profile for these models is plotted in Fig. [Ml 
which should be compared to the "natural" Fig. [T] and 
the more contrived toy model depicted in Fig. [TTJ These 
models can be constrained using the methods proposed 
in this Section. In the face of a negative result, as above, 
we'd be able to assume we are in the quasi-Newtonian 
regime, to conclude, after a similar argument, the coun- 
terpart of 



(5F. 



47r 



{5ti)¥ 



N 



K 

An 



An flo 



' N\ 



' N 



(57) 



We can now constrain parameter n as before, with the 
result plotted in Fig. [T5l As in Fig. [T3] we have plotted 
the value of n (not to be confused with the parameter 
used there) for which the SNR turns out 1, for a given 
noise and impact parameter. A negative result from LPF 
would therefore require n to be smaller than this value, 
which is therefore an upper bound. 



VII. CONCLUSIONS 

In this paper we showed how a LPF saddle flyby would 
either detect MOND to a high SNR or rule it out, if 
not comprehensively, at least to a large extent. The for- 
mer conclusion could be expected. Even though in Sec- 
tions mil and IIVI we provided rigorous and quantitative 
SNR estimates, the high levels forecast can be understood 
with a "back of the envelope" calculation. The exercise 
highlights an uncanny coincidence. The accelerometer 
aboard LPF has a non-white noise profile, dipping in the 
region of the mHz, i.e. in the rough time scale of minutes. 
The motivation for this design lies in the gravitational 
wave signals to be targeted by LISA. It just so happens 
that the MONDian bubbles of anomalous tidal stresses 
around the Earth-Sun-Moon saddles have a length scale 
of the order of a hundred kilometers. Anything free- 
falling in the Earth-Moon region has a typical speed of 
the order of Ikm/s. Thus, the time scale for crossing a 
MONDian bubble will be of the order of minutes: just 
where the instrument performance is optimal. This is 
a remarkable coincidence. Scribbling on the back of an 
envelope, using the expression for the SNR of a noise- 
matched filter and the order of magnitude of the stresses 
and noise, promptly reveals double figure SNRs. 
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FIG. 15: Contours of the power n for models (I55p required in 
order to obtain SNR= 1, for different noise levels and impact 
parameters up to b = 600 km. These are to be seen as an 
upper bound on parameter n, should a negative result from 
a LPF saddle experiment be found, for a given experimental 
set up. 



The question then arises as to how generic this conclu- 
sion is, or conversely, should a negative result be found, 
how thoroughly have we ruled out MOND. In Section HIl 
we scanned the full array of MONDian theories, show- 
ing that in the non-relativistic regime they fall into only 
3 categories (which we labelled type I, IIA/B and III). 
Even if we restrict ourselves to type I theories (a class 
containing the vast majority of these theories, including 
TeVeS), we benefit from the leverage of a free function, fi. 
Could MONDologists use /i to survive a negative result? 

In Sections |V] we examined /i functions on offer in the 
literature and laid down criteria for reasonable fi based 
on astrophysical usefulness, viability in the face of con- 
straints, and naturalness. We found that once these crite- 
ria are taken into account the size of the MOND bubble, 
rg, is fixed. There are exceptions to this rule; e.g. di- 
verging /x functions, but these may be subject to other 
constraints (they will be examined elsewhere). Predic- 
tions for what happens inside the bubble are also model 
independent; however the tidal stress anomalies outside 
the bubble depend on the transient from MONDian into 
Newtonian regime, with a fall-off which is indeed model 
dependent. Thus, for impact parameters smaller than tq 
the predicted SNRs are robust, and do not change sub- 
stantially with the model. For the currently expected b 
(around 50 km, with tq '^ 380 km) this is indeed the 
case. 

Therefore the only way for MOND to wriggle out of 
a negative LPF result would be to change the bubble 
size rg. This can only be accomplished with "designer" 
/i- functions. If fi is allowed to have two scales and two 



power-laws away from its Newtonian value of 1, then it 
is possible to bypass a negative LPF result. Even for 
undemanding noise levels and impact parameters, the 
intermediate power becomes very contrived. Similarly 
fine-tuned functions have been proposed in the literature. 
In Section IVI A[ we showed how LPF could be used to 
constrain them. The point remains that one would have 
to bend backwards to accommodate a negative result, 
within type I theories. There are exceptions to this rule 
(for example diverging /i functions) and these derserve 
further attention. 

Although we didn't present quantitative results for 
type II theories, the same conclusions apply qualitatively 
to type IIB theories (but not type II A theories). If G is 
renormalized and the free function v is chosen to produce 
the same phenomenology as type I theories (in particu- 
lar with regards to Gaen and MONDian behaviour), then 
the MOND bubble has the same size, and the anomalous 
tidal stresses are of the same order. As explained in Sec- 
tion ini in both types of theory MONDian behaviour is 
due to an extra field </>, and if one attends simultaneously 
to Giien ~ G and (f> ^ ^n for cat « oq, then MONDian 
behavior in should be triggered at the same Newtonian 
acceleration a at — a^^^ ^ oq. This implies a MONDian 
bubble of the same size tq. Furthermore the (also v- 
independent) effects inside the bubble are different from 
type I predictions, but stronger. Type II theories don't 
have a curl field (in the sense define above), a feature 
which softens the anomalous tidal stresses in type I the- 
ories [Tl|. A detailed quantitative prediction for type II 
theories is currently being investigated (see [2^; also we 
reference here a paper on the matter which has appeared 
since the present paper was submitted J31|.) 

However, it may be that the relativistic mother the- 
ory is set up in such a way that the cosmological and 
non-relativistic G coincide, the case of type IIA theories. 
In this case the MOND bubble around the saddle is very 
small. Likewise type III theories (or, rather, its single rel- 
ativistic realization, the Einstein- Aether theory) produce 
effects around saddles which are unobservable with cur- 
rent technology. In these theories G is not renormalized 
and a*"^ = ag, so that the MOND bubble is a few me- 
ters across. Remarkably, solar system tests are extremely 
constraining upon type III theories, due to the so-called 
external field effect [^. By contrast solar system effects 
for type I and IIB theories are suppressed by a factor of 
k/Att. Thus saddle tests and planetary orbits seem to be 
complementary in constraining MONDian theories. 

We close by noting that we could, of course, detach our 
considerations entirely from the MOND paradigm (as an 
alternative to dark matter) , and regard these theories for- 
mally as a class on alternative theories of gravity (see (33 
for an extensive review) . It is remarkable that only three 
classes of theories emerge in the non-relativistic regime, 
which we labelled type I, II and III in Section |lll We 
could then view k, and ao as free parameters, converting 
a LPF saddle fiyby into a constraint or a detection in 
this space. We are currently working on this alternative 



17 



approach. 
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